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ABSTRACT 

We report high resolution 12 CO(J = 2-1), 13 CO(J = 2-1), and 12 CO(J = 3-2) imaging of the 
Seyfert 1/starburst ring galaxy NGC 1097 with the Submillimeter Array to study the physical and 
kinematic properties of the 1-kpc circunmuclear starburst ring. Individual star clusters as detected 
in the HST map of Paa line emission have been used to determine the star formation rate, and 
are compared with the properties of the molecular gas. The molecular ring has been resolved into 
individual clumps at GMA-scale of 200-300 pc in all three CO lines. The intersection between the 
dust lanes and the starburst ring, which is associated with orbit-crowding region, is resolved into 
two physically/kinematically distinct features in the l'.'5xl'/0 (105x70 pc) 12 CO(J = 2-1) map. The 
clumps associated with the dust lanes have broader line width, higher surface gas density, and lower 
star formation rate, while the narrow line clumps associated with the starburst ring have opposite 
characteristics. Toomre-Q value under unity at the radius of the ring suggests that the molecular ring 
is gravitationally unstable to fragment at the scale of the GMA. The line widths and surface density 
of gas mass of the clumps show an azimuthal variation related to the large scale dynamics. The star 
formation rate, on the other hand, is not significantly affected by the dynamics, but has a correlation 
with the intensity ratio of 12 CO (J = 3-2) and 12 CO(J = 2 1), which traces the denser gas associated 
with star formation. Our resolved CO map, especially in the orbit-crowding region, for the first time 
demonstrates observationally that the physical/kinematic properties of the GMAs are affected by the 
large scale bar-potential dynamics in NGC 1097. 

Subject headings: Galaxies: individual (NGC 1097) - Galaxies: ISM - Galaxies: Seyfert - Galaxies: 
starburst - radio lines: ISM - Submillimeter: ISM 



1. INTRODUCTION 

NGC 1097 [SB(s)b; Ide Vaucouleurs et al.1 [199~l| is a 
nearby (D = 14.5 Mpc; 1" = 70 pc, |TjJlyJ [L98|| barred 
spiral galaxy. A pair of dust lanes are located at 
the leading edges of the major bar. A radio contin- 
uum image at 1.465 G Hz shows faint ridg es coincid- 
ing with the dust lanes (jHummel et al.l [1987) . The nu- 
cleus is thoug ht to be a transition object fr om LINER 
to Seyfert 1 (jStorchi-Bergmann et al.1 120031 ). Detailed 
studies on the nucleus show morphological and kine- 
matic evidences of the nuclear spirals on the order 
of 30 pc, and was inte rpreted as part of the fueling 
chain to the very center (iFathi et al.ll20ul ID avies et all 
[20091 Ivan de Ven & FatLill201flfrirGC 10 97 is also an 
IRAS bright galaxy ( [Sanders et all 120031). The con- 
tribution of large amount of IR flux a rise from its 1 
kpc-c i rcumnuclear starburst r i ng (e.g.. iHummel et al.1 
[19871: iTelesco fc Gatlevl [l98l IKotilainen et al.l 120000 . 
The starburst ring hosts "hot-spots" compo sed with su- 
per s tar clusters identified in HST images (jBarth et al.1 
1995), and was suggested to have an instantaneous 
burst of star formation which occurred ~ 6-7 Myr ago 
(jKotilainen et al.ll2000D . 

The molecular gas of NGC 1097 in the nuclear region 



has been previously mapped in the dense gas tracer of 
HCN(J = 1-0), lo w excitation lines of 12 CO(J = 1-0) 
and 12 CO (J = 2-1) (IKohno et al][2003l:IHsieh etalMM . 
hereafter Paper I). These maps show a central concen- 
tration coinci dent with the peak o f the 6-cm radio con- 
tinuum core (jHummel et al.l 119871 ) , as well as a molec- 
ular ring coincident with the starburst ring. A pair 
of molecular ridges coincident with the dust lanes are 
also detected, and show non-circular mot ions, possibly 
caused by the bar-potential dynamics ( e . g . . I Athanassoulal 
ll992allbT ). The molecular ring has a typical warm tem- 
perature (7k ~ 100 K) and denser gas (nn 2 ~ 10 3 cm" 3 ) 
consi s tent with the star burst environments (|Wild et al.l 
Il992t lAalto et al.lll995f) . The molecular ring exhibits 
a twin-peak structure in the 4"-10" resolution inter- 
fcrometric CO and HCN maps, where a pair of molec- 
ular concentrations are located in the intersection of 
the molecular dust lanes and the star forming ring. 
Its orientation is nearly perpendicular to the stellar 
bar. The twin-peak has higher H2 column density than 
the surrounding ring, and similar features have been 
seen in other barred galaxies, and can be explained 
by the crowding of gas streamlines (e.g.. iKennev et al.1 
[l99llRevnaud fe Downei[l997l: iKohho et al.lll999D . The 
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gas flow gradually changes direction and migrates to- 
ward t he center of the galaxy to accumu l ate to form 
a rin g (jSchwarzl H98l IAthanassouIalll992bl; iPiner et al.l 
I1995T ). Subsequent enhanced star formation may oc- 
cur through grav itational fragmentation stochastically 
(jElmegreenl I1994T ) . or dynamically driven by collision 
of molecular clouds (|Combes fc Gerinlll985D . or alterna- 
tively originated from the shock compressed gas near the 
conta ct point of the dust lanes and the ring (|Boker et al.l 
2008). Another intriguing topic is thus whether the oc- 
currence of the circumnuclear starburst ri ng would pro- 
hibit or boost its nuclear activities (e.g.. iTelesco et al.l 



filial IScovilie et al.lUl^lHeckmaiJl99HlHo et al.ll!997fl T 
High spatial/kinematic resolution observations of 
molecular lines are essential to study the circumnuclear 
ring structures, since they are the sites of star formation 
and respond to the large scale dynamics. NGC 1097 is 
one of the best example to study the circumnuclear ring 
for its typical structures of dust lanes, starburst ring, 
and nuclear activities. In order to study the physical 
and kinematic properties of the starburst ring, especially 
in the twin-peak region, we have now obtained higher 
resolution 12 CO(J = 2-1) (l"5x l'.'O), 13 CO(J = 2-1) 
(l"8xl"4), and 12 CO(J = 3-2) (3'.'5x2'.'l) maps down to 
100 pc. By virtue of the high angular resolution multi-J 
lines, we derive the fundamental properties of molecular 
gas as well as star formation in the ring in order to give 
a comprehensive view of this system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. SMA observations 

We observed N GC 1097 with the Submillimeter Arrajfl 
(|Ho et al.ll2004f) at the summit of Mauna Kea, Hawaii. 
The array consists of eight 6-m antennas. Four basic 
configurations of the antennas are available. With the 
compact configuration, two nights of 12 CO(J = 2-1) data 
were obtained in 2004 (Paper I). To achieve higher spa- 
tial resolution, we obtained two further nights of 12 CO(J 
= 2-1) data with the extended and the very extended 
configurations in 2005. To study the excitation of the 
gas, we also obtained one night of 12 CO (J = 3-2) data 
in 2006 with the compact configuration. All the observa- 
tions have the same phase center. We loc ated the phase 
cente r at the 6-cm peak of the nucleus ([Hummel et al.l 
1987). Detailed observational parameters, sky condi- 
tions, system performances, and calibration sources are 
summarized in Table [1] 

The SMA correlator processes two IF sidebands sep- 
arated by 10 GHz, with ^2 GHz bandwidth each. The 
upper /lower sidebands are divided into slightly overlap- 
ping 24 chunks of 104 MHz width. With the advantage of 
this wide bandwidth of the SMA correlator, the receivers 
were tuned to simultaneously detect three CO lines in the 
230 GHz band. The 12 CO (J = 2-1) line was set to be 
in the upper sideband, while the 13 CO(J = 2-1) and the 
C 18 0(J = 2-1) lines were set to be in the lower sideband. 
For the 12 CO(J = 3-2) line, we placed the redshifted fre- 
quency of 344.35 GHz in the upper sideband. 

The SMA is a joint project between the Smithsonian Astro- 
physical Observatory and the Academia Sinica Institute of Astron- 
omy and Astrophysics and is funded by the Smithsonian Institution 
and the Academia Sinica. 



We calibrated the SMA data with the MIR-IDL soft- 
ware package. The detailed calibration procedures of the 
data of compact configuration were described in Paper 
I. For the extended and very extended configurations, 
where Uranus is resolved quite severely, we observed two 
bright quasars (3C454.3 & 3C111), and adopted a similar 
bandpass/flux calibration method as described in Paper 
I. At 345 GHz, we used both Uranus and Neptune for 
the flux and bandpass calibrators. 

Mapping and analysis were done with the MIRIAD 
and NRAO AIPS packages. The visibility data were 
CLEANed in AIPS by task IMAGR. We performed the 
CLEAN process to deconvolve the dirty image to clean 
image. The deconvoluti on procedure w as adopted with 
the Hogbo m algorithm (|Hogboml I1974D and the Clark 
algorithm (jClarkl [1980) . We used a loop gain of 10% 
and did restrict the CLEAN area by iterative examin- 
ings. The CLEAN iterations were typically stopped at 
the 1.5cr residual levels. The number of CLEAN com- 
ponents in individual channel maps was typically about 
300. 

All of the 230 GHz visibility data were combined to 
achieve a better uv coverage and sensitivity. Our 230 
GHz continuum data were constructed by averaging the 
line-free channels. Due to the sideband leakage at 230 
GHz, a limited bandwidth of 1.3 GHz and 0.5 GHz were 
obtained respectively in the upper and lower sidebands 
to make a continuum image. Weak continuum emission 
at 230 GHz with a peak intensity of 10 mJy beam -1 
was detected at about 4cr at the southern part of the 
ring. The la noise level of the continuum emission at 
345 GHz averaged over the line-free bandwidth of 0.4 
GHz is 6 mJy beam -1 , and there is an ~ 3er detection 
in the nucleus and in the ring. In this paper we did not 
subtract the continuum emission since it is too faint as 
compared to the noise level in the line maps. The spec- 
tral line data of 12 CO(J = 2-1) and 12 CO(J = 3-2) were 
binned to 10 km s -1 resolution. As the 13 CO(J = 2-1) 
emission is fainter, those data were binned to 30 km s -1 
resolution to increase the S/N. The C 18 0(J = 2-1) line 
emission was detected in the lower sideband of the 230 
GHz data. However, the leakage from the 12 CO (J = 2- 
1) line in the upper sideband was significant. Therefore, 
in this paper we will not make use of the C ls O(J = 2- 
1) data for further analysis. In this paper, we present 
the maps with natural weighting. The angular resolu- 
tion and rms noise level per channel are l'/5 x l'.'O and 26 
mJy beam" 1 (400 mK) for the 12 CO(J = 2-1), 1'.'8 x 1'.'4 
and 30 mJy beam" 1 (320 mK) for the 13 CO(J = 2-1) 
data, and 3'.'5 x 2'.'1 and 35 mJy beam -1 (47 mK) for the 
12 CO(J = 3-2) data. We used AIPS task MOMNT to con- 
struct the integrated intensity-weighted maps. The task 
would reject pixels lower than the threshold intensity, set 
to be 2.5-3<7 in the CLEANed cube after smoothing in 
velocity and spatial direction. The smoothing kernels are 
30 km s _1 in the velocity direction, and a factor of two 
of the synthesized beam in the spatial direction in the 
maps we made. 

In the following line ratio analysis, we use the maps 
truncated to the same uv coverage as in Sect. 13.4.11 



2.2. HST NICMOS Pan image 



CO(J = 2-1) Imaging of Seyfert 1 Galaxy NGC 1097 
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As one of the sample galaxies in the HST NICMOS 
survey of nearby galaxies (PI: Daniela Calzetti, GO: 
11080), NGC 1097 was observed on November 15, 2007 
by HST equipped with NIC3 camera. NIC3 images 
have a 51"x51" field of view and a plate scale of 0'/2 
with an undersampled PSF. In this survey, each observa- 
tion consists of images taken in two narrowband filters: 
one centered on the Paa recombination line of hydro- 
gen (1.87 /Ltm) (F187N), and the other on the adjacent 
narrow-band continuum exposure (F190N), which pro- 
vides a reliable continuum subtraction. Each set of ex- 
posures was made with a 7-position small (< 1" step) 
dithered. Exposures of 160 and 192 seconds per dither 
position in F187N and F190N, respectively, reach a la 
detection limit of 2.3 xl0~ 16 erg s _1 cm~ 2 arcsec~ 2 in 
the continuum-subtracted Paa image of NGC 1097. 

Using the STSDAS package of IRAF, we removed the 
NICMOS Pedestal effect, masked out bad pixels and cos- 
mic rays, and drizzled the dither images onto a finer (O'.'l 
per pixel) grid frame. The resultant drizzled F190N im- 
age is then scaled and subtracted from its F187N peer af- 
ter carefully aligning to the latter using foreground stars. 
The residual shading effect in the pure Paa image is re- 
moved by subtracting the median of each column. This 
strategy works well for the NGC 1097 data which con- 
tain relatively sparse emission features. The PSF of the 
final Paa image has a 0'.'26 FWHM. Further details of 
th e data reduction and image processing are described 
in lLiu et al.l (|2010aj ). 

3. RESULTS 
3.1. Morphologies of the molecular gas 

The integrated intensity maps of 12 CO(J = 2-1), 
12 CO(J = 3-2) and 13 CO(J = 2-1) lines are shown in 
Figure [T] The maps show a central concentration and 
a ring-like structure with a radius of 700 pc (~10"). 
The gas distribution of 12 CO (J = 3-2) map is similar 
to that of the 12 CO(J = 2-1) map, where the central 
concentration has a higher integrated intensity than the 
ring. The 13 CO(J = 2-1) map, on the other hand, shows 
comparable intensity between the ring and the nucleus. 
Comparing with the previous observations (Paper I) , the 
molecular ring and the central concentration have been 
resolved into individual clumps, especially for the twin- 
peak structure in the molecular ring, with the higher 
resolution 12 CO(J = 2-1) map. We also show the peak 
brightness temperature map of the 12 CO(J = 2-1) emis- 
sion in Figure [T] made by the AIPS task SQASH to ex- 
tract the maximum intensity along the velocity direction 
at each pixel. Note that the ring and the nucleus have 
comparable brightness temperatures. 

The position of th e AGN (6-cm radio continuum core; 
iHummel et al.llT987[ ). which is assumed to be the dynam- 
ical center of the galaxy, seems to be offset by 0'.'7 north- 
west to the central peak of the integrated CO maps. This 
is also seen in our previous results and we interpreted it 
as a result of the intensity weighting in the integration 
(Paper I). To confirm if the dynamical center as derived 
from the 12 CO (J = 1-0) emission (|Kohno et al.ll200l is 
consistent with that of 12 CO(J = 2-1), we will have some 
more analysis in Sect. 13.51 

In Figures H and H we show the 12 CO(J = 2-1) and 
12 CO (J = 3-2) channel maps overlaid on the archival 



HST I-band image (F814W) to show the gas distribu- 
tion at different velocities. We show these images at the 
resolution of 40 km s , but we use the resolution of 10 
km s _1 for actual analysis. The astrometry of the HST 

1- band image was corrected by the known positions of 
the foreground 19 stars in the catalog of USNO-A 2.0. 
After astrometry correction, we found the intensity peak 
of the nucleus of HST image has a 0'/8 offset northeast to 
the position of AGN, this may be due to extinction, or 
inaccurate astrometry. We can not rule out these factors. 
However this offset is within the uncertainty of the CO 
synthesized beam. Both 12 CO transitions show emission 
with a total velocity extent of ~550-600 km s _1 at the 2a 
intensity level. The western molecular ridge/arms (coin- 
cident with the dust lane in the optical image) joins the 
southwestern part of the ring from -162 km s _1 to 37 
km s _1 , and the eastern molecular ridge joins the north- 
east ring from -42 km s~ 1 to 157 km s^ 1 . We find that 
the velocity extent of these ridges is ^200 km s^ 1 at the 
2a intensity level. In between the molecular ring and 
the nuclear disk, there are also extended ridge emission 
connecting the nuclear disk and the ring in both lines. 
This ridge emission is more significant in the 12 CO(J = 

2- 1) than in the 12 CO(J = 3-2) map. However, this 
may be an angular resolution effect. The ridge emis- 
sion looks more significant in the 12 CO(J = 2-1) map, 
but the ridge emission does not have enough flux density 
(Jy/beam) to reach the same S/N in the high resolution 
image, which means it has extended structure, i.e., the 
faint ridge emission in the high resolution map is more 
or less resolved. 

In Figure [U we show the 12 CO(J = 2-1) integrated 
intensity map overlaid on the HST I-band (F814W) 
archival image to compare the optical and radio mor- 
phologies. The Paa image (near-IR) is also shown in 
Figure 2] The optical image shows a pair of spiral arms 
in the central 1 kpc region. The arms consist of dusty 
filaments and star clusters. Two dust lanes at the lead- 
ing side of the major bar connect to the stellar arms. 
Dusty filaments can also be seen to be filling the area 
between the spiral arms and the nucleus. The 12 CO (J 
= 2-1) map (l'.'5xl'/0) shows a central concentration, 
a molecular ring and molecular ridges with good gen- 
eral correspondence to the optical image. However, there 
are significant differences in a detailed comparison. Al- 
though the starburst ring looks like spiral arms in the 
optical image, the molecular gas shows a more complete 
ring, and the molecular ridges are seen to join the ring 
smoothly. Inside of the stellar ring, the major dust lanes 
are offset from the molecular ring. Assuming that the 
12 CO emission is faithfully tracing the total mass, the 
prominent dark dust lanes are not significant features 
while they join the molecular ridges as the edge of the 
ring. The stellar ring or star formation activities are then 
correlated with the peaks in the total mass distribution. 
The central molecular concentration corresponds in gen- 
eral with the stellar nucleus, but is offset to the south 
of the stellar light. The nuclear molecular distribution is 
quite asymmetric with lots of protrusions at the lower in- 
tensity contours, which may correspond to gas and dust 
filaments which connect the nuclear concentration to the 
molecular ring. 
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3.2. Properties of the molecular gas 

The molecular ring has been resolved into a complex 
structure of compact sources immersed in a diffuse emis- 
sion with lower surface brightness as shown in Figure [TJ 
At our resolution (100 pc), the contrast between clumps 
and diffuse emission is not high. Hence, the clumps are 
possibly connected to each other via the diffuse emis- 
sion. With such a physical configuration, it is difficult 
to uniquely isolate the individual clumps by some clump 
finding algorithm fe.g. JWilliams. et al.lll994T ). Moreover, 
the detection of clumps is dependent on the angular res- 
olution which is available. In this paper, we select the 
peaks of the main structures in the 12 CO(J = 2-1) in- 
tegrated intensity map, in order to locate the individual 
clumps within the ring and the molecular ridges. We de- 
fine the clumps located outside of the 10" radius as the 
dust lane clumps, and those within the 10" radius as the 
ring clumps. 

The typical size of Giant Molecular Cl ouds (GMCs) 
is on the order of a few to few tens pc (jScoville et al.l 
[19871 ). The size of the clumps in the ring of NGC 1097, 
as detected with our synthesized beam, is at least —200 
pc. We are therefore detecting molecular clumps larger 
than GMCs, most likely a group of GMCs, namely Gi- 
ant M olecular Cloud Associations (GMAs; Vog el et"aH 
Il998h . Here we still use the term "clumps" to describe 
the individual peaks at the scale of GMA. We expect 
that the clumps would be resolved into individual GMCs 
when higher angular resolution is available. In Table [H 
we show the quantities measured within one synthesized 
beam to study the kinematic properties with a high res- 
olution (i.e., one synthesize beam) in the following sec- 
tions. The observed peak brightness temperatures of in- 
dividual clumps are in the range from —2 - 8 K. 

To show the general properties of the GMAs, we also 
measured the physical properties integrated over their 
size in Table [3] We define the area of the GMA by mea- 
suring the number of pixels above the threshold intensity 
of 5(7 in the 12 CO(J = 2-1) integrated intensity map. 
We calculated the "equivalent radius" if the measured 
areas are modeled as spherical clumps. The results are 
reported in Table [3] together with the resulting Mh 2 in- 
tegrated over the area. The derived values of A/h 2 are 
therefore larger than that in Table O which only mea- 
sured the mass within one synthesized beam at the in- 
tensity peak. The method to derive the Mh 2 will be de- 
scribed in Sect. 13.2^2"! Several factors are essential to con- 
sider for fair comparisons of our GMAs with other galax- 
ies, such as beam size, filling factor, etc. A rough com- 
parisons show that our GMAs in the starburst ring have 
physical extent of —200-300 pc, which have similar order 
with the GMAs in other galaxies fe.g.lRand &; Kulkarnil 
[19901 iTosaki et all 120071 iMuraoka et al.ll2009D . 

3.2.1. Spectra of the molecular clumps 

We show in Figure \5\ the spectra of the individual 
clumps selected above at their peak positions (i.e., within 
a synthesized beam). Most of the clumps show single 
gaussian profiles. The gaussian FWHM line widths of 
the clumps are determined by gaussian fitting, and are 
listed in Tabled 

The line widths we observed are quadratic sums of the 
intrinsic line widths and the galactic motion across the 



synthesized beam, and the relation can be expressed as 

^obs = a rot + <J i 2 nt> (!) 

where er b s is the observed FWHM line width, er rot is the 
line width due to the galactic rotation and the radial mo- 
tions within the beam, and the a- lri t is the intrinsic line 
width. The systematic galactic motion would be negligi- 
ble if the spatial resolution is small enough. We subtract 
cr ro t as derived from a circular rotation model by fitting 
the rotation curve in Sect. 13.51 The derived intrinsic line 
width of the clumps are reported in Tabled) In general, 
there is —10-20% difference between the observed and 
intrinsic line widths, which suggests that the galactic cir- 
cular motions do not dominate the line broadening at this 
high resolution. However, the large "intrinsic" velocity 
dispersion thus derived by subtracting circular motions 
could still be dominated by non-circular motions, espe- 
cially in the twin-peak region and in the dust lanes. We 
will mention this effect in the discussion. 

In Table we define the clumps based on their loca- 
tion and their velocity dispersions. The clumps in the 
ring are further designated by their velocity dispersion 
being broader or narrower than 30 km s^ 1 , and named 
respectively as Bl,..., B3, and Nl,..., Nil, respectively. 
The clumps located at the dust lanes are named Dl,..., 
D5. 

3.2.2. Mass of the clumps 

Radio interferometers have a discrete sampling of the 
uv coverage limited by both the shortest and longest 
baselines. Here, we estimate the effects due to the miss- 
ing information. In Figure [6j we convolved our newly 
combined SMA data to match the beam size of the JCMT 
data (21"), and overlaid the spectra to compare the flux. 
The integrated 12 CO(J = 2-1) fluxes of the JCMT and 
SMA data are -120 K km s^ 1 and -71 K km s -1 , re- 
spectively. The integrated 13 CO(J = 2-1) fluxes of the 
JCMT and SMA data are -15 K km s" 1 and -7 K 
km s -1 respectively. Therefore our SMA data recover 
- 60% and - 47% of the 12 CO(J = 2-1) and 13 CO(J 
= 2-1) fluxes measured by the JCMT, respectively. If 
the missing flux is attributed to the extended emission, 
then the derived fluxes for the compact clumps will re- 
main reliable. However, the spectra of the SMA data 
seem to have similar line profiles as that of the JCMT. 
Part of the inconsistency could then possibly be due to 
the uncertainty of the flux calibration, which are 15% 
and 20% for the JCMT and SMA, respectively. We note 
therefore that the following quantities measured from the 
flux will have uncertainties of at least 20% from the flux 
calibration. 

There are several ways to calculate the molecular gas 
mass. The conventional H2 / Iqo (Ago ) conversion factor 
is one of the methods to derive the gas mass assuming 
that the molecular clouds are virialized. First, we use 
the 12 CO emission to calculate the molecular H2 column 
density, Ah 2 , as 

A H2 = A co J T h dv = X co I C o, (2) 

where Ah 2 for each clump is calculated by adopting the 
Aqq conversion fact or of 3xl0 20 cm~ 2 (K km s -1 ) -1 
IjSolomon et al.|[l987[ ), Tb is the brightness temperature, 



CO(J = 2-1) Imaging of Seyfert 1 Galaxy NGC 1097 



5 



and dv is 
Table H Note that t 
4xl0 20 cm" 2 (K km 



the line width. The values are listed in 



s- 1 )- 1 (Young & Scovillc 1991; 


iDame et all 2001 




Draine et al.l 


20011: IMeier et al. 




2008). We 



adopted 3xl0 20 cm -2 (K km s -1 ) -1 to be consistent 
with Paper I. The surface molecular gas mass density 
Sh 2 and the molecular gas mass Mh 2 are thus calculated 

as, 



S H2 = 



M H2 = S H2 dO, 



(3) 



(4) 



where mn 2 and dil are the mass of a hydrogen molecule 
and the solid angle of the integrated area, respectively. 
Since this conversion factor is used for the 12 CO(J = 1 
0) emission, we assume the simplest case where the ratio 
of 12 CO(J = 2-l)/ 12 CO(J = 1-0) is unity. The H 2 mass 
of the clumps measured within one synthesized beam are 
listed in Table [2] We would get the molecular gas mass 
(M gas ) by multiplying the H 2 mass by the mean atomic 
weight of 1.36 of the He correction. 

The 12 CO(J = 2-1) is usually optically thick and only 
traces the surface properties of the clouds. Optically 
thinner 13 CO(J = 2-1) would be a better estimator of 
the total column density. We therefore calculate the av- 
erage Mjj 2 in the nucleus and the ring to compare the 
M H2 derived from both 12 CO(J = 2-1) and 13 CO(J = 
2-1) lines. In the case of 12 CO(J = 2-1), the average 
Mh 2 of the ring is about 1.7xl0 7 M within one synthe- 
sized beam, and the corresponding value for the nucleus 
is 3.4 xlO 7 Mq. However, in Paper I we derived the in- 
tensity ratio of 12 CO(J = 2-l)/ 12 CO(J = 1-0) is -2 
for the nucleus in the lower resolution map. Therefore 
the Mh 2 of the nucleus is possibly smaller than that de- 
rived above. The conventional Galactic Xco of 3xl0 20 
cm~ 2 (K km s -1 ) -1 is often suggested to be overesti- 
mated in the galactic center an d the starburst environ- 
ment by a factor of 2 to 5 (e.g., Maloncv & Black3 fi~988l : 
IMeier k, Turnerll2001| ). Therefore, our estimated M Ha in 
the ring and the nucleus might be smaller at least by a 
factor of 2. 

Assuming the 13 CO (J = 2-1) emission is optically thin, 
and t he 13 CO/H 2 abundance of lxl0~ 6 (Sol omon et al.l 
1979). We calculate the Mh 2 with excitation temper- 
ature (T ex ) of 20 K and 50 K in LTE condition. The 
Mh 2 of the nucleus averaged over one synthesized beam 
are (7.0±2.5) x 10 6 M Q and (1.3±0.4) x 10 7 M for 20 K 
and 50 K, respectively. The Mh 2 of the ring averaged 
over one synthesized beam are (5.8±2.5) x 10 6 M Q and 
(1.1±0.4) x 10 7 M Q for 20 K and 50 K, respectively. With 
the assumption of constant 13 CO/H2 abundance, if the 
T ex is < 20 K, then the conversion factor we adopted for 
the 12 CO line is overestimated by a factor of ~3 to 5. 
The overestimation is smaller (a factor of 2 to 3) if the 
gas is as warm as 50 K. We measured the 12 CO(J = 2-1) 
flux higher than 3tr to derive the total H 2 mass of the 
nucleus and the ring. The total flux of the nucleus and 
the ring are 490.2±59.3 Jy km s" 1 and 2902.2±476.8 Jy 
km s — 1 , respectively. If we adopted the intensity ratios 
of 12 CO(J = 2-1) and 12 CO(J = 1-0) are 1.9±0.2 and 



1.3±0.2 for the nucleus and the ring (Paper I), the Mh 2 
of the nucleus and the ring are (1.6±0.3) x 10 8 M© and 
(1.4±0.3)xl0 9 M Q , respectively. Thus the nucleus and 
the ring account for 10% and 90% of the Mh 2 within the 
2 kpc circumnuclear region, respectively. 
We calculate the virial mass of individual clumps by 



„ 2 < 

M v ir — — . 



(5) 



where r is the radius of the clump derived in 
Sect. I3.2[ and <7 rms is the three dimensional intrinsic 
root mean square velocity dispersion, which is equal to 
y/ (3 /8in2)cri n t ■ Since we observe the a m t in one dimen- 
sion, we need to multiply the observed af nt by 3. The 
radius is adopted with the size of the clumps assuming 
that (Ti n t is isotropic. The results are shown in Table |3l 
Note that the virial mass is for GMAs not GMCs, and we 
assume that the GMAs are bounded structures. The ra- 
tio of the virial mass M v ; r and M gas are shown in Table |3l 
We found that the narrow line clumps have M v i r /M gas 
that are more or less about unity, but are larger than 
unity for the broad line and dust lane clumps. 

We plot the general properties of the molecular gas 
mass of individual molecular clumps in Figure [7l In Fig- 
ure [7^, the histogram of the gas mass integrated over 
their size seems to be a power law with a sharp drop at 
the low mass end. This is due to the sensitivity limit, 
since the corresponding 3er mass limit is - 27 x 10 6 M Q 
for a clump with a diameter of 3'.'3 (the average of the 
clumps in Table [3]). In Figure [7Jd, the FWHM intrinsic 
line widths have a weak correlation with gas mass in the 
ring. In Figure [8^, we show the azimuthal variation of 
Sh 2 calculated at the emission peaks of the clumps in- 
tegrated over one synthesized beam. If we assume that 
NGC 1097 has a trailing spiral, then the direction of rota- 
tion is clockwise from east (0°) to north (90°). The Eh 2 
in the orbit crowding region is roughly from 0° to 45°, 
and from 180° to 225°. Note that the dust lane clumps 
typically have £h 2 similar to the narrow line ring clumps, 
which are lower than the broad line ring clumps in the 
orbit crowding region. The average Sh 2 of the narrow 
line ring and dust lane clumps is ~ 1800 M© pc~ 2 , and 
~ 2300 M© pc -2 for the broad line ring clumps. In Fig- 
ure [8p, the velocity dispersion of clumps located at the 
orbit crowding region and the dust lane are larger than 
that of the narrow line ring clumps. The average velocity 
dispersion of the narrow line clumps is ~ 50 km s _1 , and 
~ 90 km s _1 for the broad line ring/dust lane clumps. 
Given the similar line brightnesses (Table[2]) of the peaks, 
the increased £h 2 are probably the results of increased 
line widths as it is indicated in equation (2). On the 
other hand, since the Nn 2 is proportional to the number 
density of gas and line-of-sight path (i.e., optical depth), 
the higher Eh 2 of the broad line clumps may be due to 
either a larger number density, or a larger line-of-sight 
path. We will discuss these effects in the discussion. 

3.3. Young star clusters 

To check how the star formation properties are asso- 
ciated with the molecular clumps in the ring, we com- 
pare the 6-cm radio continuum, V-band, and Paa im- 
ages with the molecular gas images. The 6-cm radio 
continuum sources are selected from the intensity peaks 
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(|Hummel et all 11987ft at an angular resolution of 2'.' 5. 
The V-band clusters (< 13 mag) are s elected from the 
HST F555W image (|Barth et aliU99l . These V-band 
selected clusters have typical size of 2 pc (0V03), and ar e 
suspected to be super star clusters by Barth ct ah] (|1995D . 
Paa clusters are identified in our HST F187N image de- 
scribed in the Sect. [5] The Paa clusters are identified by 
SExtractor with main parameters of detection threshold 
of 15(J (DETECT_THRESH) and minimum number of pix- 
els above threshold of 1 (DETECT_MINAREA) . The number 
of deblending sub-thresholds is 50 for DEBLENDJJTHRESH 
and 0.0005 for DEBLEND_MINCONT. The parameters were 
chosen by a wide range of verifying and visual inspection. 

We used 6-cm radio continuum and V-band selected 
clusters for the phenomenological comparison with the 
molecular clumps. To get a high resolution and uniform 
sample of star formation rate (SFR), we use the Paa 
clusters in the following analysis. 

In the 12 CO(J = 2-1) integrated intensity map (Fig- 
ure [T]), the star clusters and radio continuum sources are 
located in the vicinity of molecular clumps within the 
synthesized beam. The distribution of the massive star 
clusters is uniform in the ring instead of showing cluster- 
ing in certain regions. The star clusters do not coincide 
with most of the CO peaks. The spatial correlation seems 
to be better in the peak brightness temperature map in 
Figure [TJ Furthermore, there are no detected star clus- 
ters and radio continuum sources in the dust lane clumps, 
namely, clumps Dl, D2, D3, and D4. 

We corrected the extinction of the Paa emission by the 
intensity ratio of Ha (CTIO 1.5 m archived image) and 
Paa. The PSFs of Ha and Paa are ~l'/0 and ~tf'3, re- 
spectively, and an additional convolving Gaussian kernel 
has been applied to both Ha and Paa images to match 
the CO beam size. The observed intensity ratio of Ha 
and Paa, (lHa/Ipaa)o, and the predicted intensity ratio 
(lHa/Ipaa)i are multiplied by the extinction as follows: 

( IE?-} = ( x lO-^^-VH'^-Kpac,^ ^ 

where -E(B-V) is the color excess, kh« and Kp aQ are the 
extinction coefficients at the wavelength of Ha and Paa, 
respectively. 

The predicted intensity ratio of 8.6 is derived in the 
case B recombination with temperature and e lectron den- 
sity o f 10000 K and 10 4 cm~ 3 , respectively (|Osterbrockl 
1989). The extinction coefficient of Ha and Paa are 
adop ted from Cardelli's extinction curve (Cardclli ct al. 
1989), with «Ha = 2.535 and Kp aQ = 0.455, respectively. 
We derive the color excess .E(B-V) using equation (6), 
and derive the extinction of Paa using the following equa- 
tion, 

A\ = k\E(B — V), (7) 

where A\ is the extinction at wavelength A. We show the 
derived values for Ap aQ and Ay in Table 01 where ny is 
3.1 at V-band. The average extinction of the clumps at 
the wavelength of 1.88 /jm (redshifted Paa wavelength) 
is about 0.6 mag, which is quite transparent as com- 
pared with the extinction at the Ha line of ~ 4 mag. We 
calculate S FRs using the Paa l uminosity based on the 
equation in lCalzetti et al.l {2008) of 

SFR(M yr- 1 ) = 4.2 x 10- 41 L PaQ (erg s" 1 ). (8) 



The SFR surface density (Esfr) is thus calculated 
within the size of the CO synthesized beam (l'.^xl'.'O). 
Note that SFR cannot be determined in the dust lane 
clumps since we do not detect significant star clusters in 
both Paa and Ha images. The low Esfr m the broad 
line ring and dust lane clumps seems to be unlikely due 
to the extinction, since the extinction of the broad line 
ring clumps are similar to that of the narrow line ring 
clumps, and the Eh 2 of the dust lane clumps are simi- 
lar to that of the narrow line ring clumps. To compare 
the star formation activities with the properties of the 
molecular gas, we measured the Esfr at the position 
of each clump (Table 2]). We show the correlation of 
Esfr an d the Eh 2 of the molecular clumps in Figure [9^i. 
This plot shows very little correlation. In Figure [TOl 
we overlay our data on the plot of Esfr and Eh 2 used 
in Kennicutt (1998) to compare the small and the large 
scale star formation. The average number follows the 
Kennicutt-Schmidt correlation closely. However, we have 
eit her lower Esfr. o r higher Eh 2 than the global values 
in iKennicuttl (|l998). This might be because our spatial 
resolution is smaller than for their data. Recent inves- 
tigations have shown that the power scaling relationship 
of the spatially-resolved S chmidt-Kennicutt law remains 
valid in t he sub-kpc scale (|Bigiel et al.ll2008l), to ~200 pc 
in M51 dLiu et al.l lMObl) and M33 dVerlev et alJ l2010t 
IBigiel et all 12010ft. but becomes invalid at the scale of 
GMC/GMAs ()Onodera et al.ll2010T ) because the scaling 
is overcome by the large scatter. The absence of a corre- 
lation in our 100 pc study is thus not surprising because 
even if the Schmidt-Kennicutt law is s till valid, the scat - 
ter is expected to as large as ~0.7 dex (|Liu et a l. 2010b), 
larger than the dynamical range of the gas density in 
Figure Other possibilities to explain the inconsis- 
tence are the uncertain conversion factor mentioned in 
Sect. 13.2.21 The conversion factor is likely to be overesti- 
mated in the galactic center and starburst region. More- 
over, regarding to the Schmidt-Kennicutt law was de- 
rived from the global galaxies that might be dominated 
by disk GMCs, our nuclear ring might not follow the 
same relation for its particular physical conditions. 

As for the distribution of Esfr in the ring, we find that 
Esfr is low in the broad line clumps in Figure ^jp, but 
do not have an obvious systematic azimuthal variation 
as Eh 2 or intrinsic line width in Figure [H In general, 
Esfr is higher in the northern ring than in the southern 
ring. This distribution, as an average quantity, shows no 
strong dependence on location within the ring. 

3.4. Physical conditions 
3.4.1. Intensity ratio of multi-J CO lines 

We compare the intensity ratio of the different CO lines 
on the same spatial scales by restricting the data to the 
same iw-range from 7.3 kA to 79.6 kA for the 12 CO(J 
= 2-1), 13 CO(J = 2-1) and 12 CO(J = 3-2) lines. The 
matched beam size of all maps is 3'.'25x2'.'55. We cor- 
rected for the primary beam attenuation in the maps. 
We measured the line intensities of individual clumps 
in the uw-matched integrated intensity maps, and calcu- 
lated the intensity ratios in Table [5] The uu-matched low 
resolution maps do have some beam smearing effects on 
the spectra. However, an examination of the line profiles 
and attempts to correct for line smearing did not affect 
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the derived line ratios to within the experimental errors. 

We estimate the densit y and temperature of the c lumps 
with the LVG analysis (jGoldreich fc Kwanl I1974D in a 
one-zone model. The collision rates of CO are from 
iFlower fc Laimavl (I1985D f o r tem peratures from 10 to 250 
K. and from IMckee et all (fl982T) for 500 to 2000 K. We 
assume 12 CO and 13 CO abundances with respect to H 2 
of 5xl0 -5 and lxlO -6 , with the observed velocity gra- 
dient of ~ 1 km s _1 pc -1 of the ring. We determined the 
velocity gradient in the Paper I by the PV diagram, and 
it is consistent in this paper. The average ratio of the 
narrow and broad line clumps are used. Clumps N4, Bl, 
Dl, D2, D4, D5 are excluded in the average ratio because 
of their large uncertainty in R43 . Therefore the average 
R32 and R13 of the narrow line clumps are 1.00±0.02 and 
9.90±2.11, respectively. The average R32 and R13 of the 
broad line clumps are 0.72±0.01 and 9.55±1.56, respec- 
tively. With the constraint of the intensity ratios within 
the uncertainty, the estimated temperature and density 
of the narrow line clumps are >250 K and (4.5±3.5) x 10 3 
cm~ 3 . The broad line clumps have temperatures of 
45 ± 15 K and density of (8.5 ± 1.5) x 10 2 cur 3 . The 
predicted brightness temperature (Tb) is ~100 K for the 
narrow line clumps and ~20 K for the broad line clumps. 
However, it seems to be inconsistent with the high/low 
Sh 2 an d low/high number density in the broad/narrow 
line clumps if we assume a constant scale height for the 
clumps. The solution may be a smaller beam filling fac- 
tor for the narrow line clumps. 

In Figure the R32 values have a positive correlation 
with Ssfr- In Figure IHJi, similar to Ssfr, R32 is slightly 
lower in the broad line ring clumps and does not show 
any systematic pattern in the azimuthal direction. 

3.5. Kinematics 

Figure ITTa is the intensity weighted isovelocity map of 
12 CO(J = 2-1). The gas motion in the ring appears to be 
dominated by circular motion, while it shows clear non- 
circular motions in the 12 CO(J = 1-0) map as indicated 
by the S-shape nearly parallel to the dust lanes. As we 
discussed in Paper I, the non-significant non-circular mo- 
tion in the 12 CO(J = 2-1) maps is perhaps because the 
dust lanes are not as strongly detected in 12 CO (J = 2-1) 
line, along with the fact that they are closer to the edge 
of our primary beam, or the non-circular motion is not 
prominent at the high spatial resolution. The circumnu- 
clear gas is in general in solid body rotation. The velocity 
gradient of the blueshifted part is slightly steeper than 
the redshifted part. We also show the intensity weighted 
velocity dispersion map in Figure [TTb . As we mentioned 
above, the velocity dispersion is larger in the twin-peak 
region, and lower in the region away from the twin-peak 
region. 

The dynamical c enter of NGC 1097 was derived by 
IKohno et all ()2003D in their low resolution 12 CO (J = 1 
0) map. With our high resolution 12 CO(J = 2-1) map, 
we expect to determine the dynamical center more accu- 
rately. We use the AIPS task GAL to determine the dy- 
namical center. In the task GAL, 12 CO(J = 2-1) intensity- 
weighted velocity map (Figure[JTJ is used to fit a rotation 
curve. The deduced kinematic parameters are summa- 
rized in Table 6. We use an exponential curve to fit the 
area within 7" in radius. The observed rotation curve 



and the fitted model curve are shown in Figure [T^J From 
the fitted parameters, we find that the offset (~ 0'.'3) of 
the dynamical center with respect to the position of the 
AGN is still within a fraction of the synthesized beam 
size. The derived V sys has a difference of ~5 km s^ 1 be- 
tween 12 CO(J = 1-0) and 12 CO(J = 2-1) data, which is 
less than the velocity resolution of the data. Upon ex- 
amining the channel maps of the 12 CO(J = 2-1) data, 
we find that the peak of the nuclear emission is almost 
coincident with the position of the AGN with an offset 
of 0'.'3. We therefore conclude that the position offset in 
the integrated intensity map, as mentioned in Sect. 13.11 
is due to the asymmetric intensity distribution. 

4. DISCUSSIONS 
4.1. Molecular ring of NGC 1097 
4.1.1. Twin-Peak structure 



In the low resolution CO maps (Paper I, IKohno et al.l 
120031) . NGC 1097 shows bright CO twin-peak structure 
arising at the intersection of the starburst ring and the 
dust lanes. The >300 pc resolution CO data show that 
the barred galaxies usually have a large amount of centra l 
concentration of mo lecular gas ([Sakamoto et al.lll999l ). 
iKennev et al.l (|1992l) found that in several barred galax- 
ies which host circumnuclear rings (M101, NGC 3351, 
NGC 6951), the central concentrations of molecular gas 
were resolved into twin-peak structures when resolution 
of ~200 pc is attained. A pair of CO intensity concen- 
trations are found in these cases, in the circumnuclear 
ring, at the intersection of the ring and the dust lane. 
Their orientation is almost perpendicular to the major 
stellar bar. The twin-peak structure can be attributed 
to the orbit crowding of inflowing gas stream lines. The 
gas flow changes from its original orbit (the so called x\ 
orbit) when it encounters the shocks, which results in 
a large deflection angle and migrate to new orbit (the 
so called 2; 2 orbit). The gas then accumulates in the 
family of the cc 2 orbits in the shape of a ring o r nuclear 
spirals (|Athanassoulalll992aHPiner et al.iri995| ). Intense 
massive star formation would follow in the ring/nuclear 
spiral once the g as becomes dense enough to collapse 
(|Elmegreenlll994h . 

In our 100 pc resolution CO map, the starburst 
molecular ring is resolved into individual GMAs. In 
the orbit crowding region, we resolve the twin-peak 
into broad line clumps associated with the curved dust 
lanes. The narrow line clumps are located away from 
the twin-peak and are associated with star forma- 
tion. This kind of "spectroscopic components" were 
also shown in several twin-peak galaxies at the intersec- 
tion o f dust lanes and circum nuclear ring, such as NGC 
1365 (jSakamoto et all I2001 . NGC 4151 ([Dumas et al.l 
I20T0 V NGC 6946 (ISchinnerer et al.l [20071), and NGC 
6951 (jKohno et al.lll999|) . However, most of the spectra 
at these intersections show blended narrow/broad line 
components, which is perhaps due to insufficient angular 
resolution. Our observations for the first time spatially 
resolved these two components toward the twin-peak re- 
gion of NGC 1097. 

It is interesting to note that the circumnuclear ring 
is nearly circular at ~42° inclination, which indicates 
its intrinsic elliptical shape in the galactic plane. The 
schematic sketch is shown in Figure [TJJ The loci of 
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dust lanes are invoked to trace the galactic shock wave, 
and their shapes are dependent on the parameters of the 
barred potential. In the case of NGC 1097 , the observed 
dust la nes resemble the theoretical studies ([Athanassoulal 
I1992H ). with a pair of straight lanes that slightly curve 
inwards in the inner ends. These findings are consistent 
with the predicted morphology from bar-driven nuclear 
inflow. The physical properties of these clumps are dis- 
cussed in the following subsection. 

4.1.2. The nature of the molecular clumps in the ring 

In Figure [T] and Table [21 the peak brightness tempera- 
tures of individual clumps are from 2 to 8 K. These values 
are lower than the typical temperatures of molecular gas 
as expected in the environment of a starburst (~ 100 K), 
and as estimated by our LVG results in Paper I and this 
work. This lower brightness temperature may be due to 
the small beam filling factor: 



/b = 



(9) 



(10) 



where /b is the beam filling factor, 6 S and #b are the 
source size and the beam size, respectively. We assume 
#s "C Oh- 7b is the brightness temperature and T c is 
the actual temperature of the clouds. We assume here 
that the source size is compact and much smaller than 
the beam. If we assume the LVG predicted Tb ~100 K 
Sect. 13.4.11 for narrow line ring clumps, and an average 
observed Tb ~5 K , then /b ~ 0.05. This suggests 9 S 
is ~20 pc or association of much smaller clumps for the 
narrow line ring clumps. In the case of the broad line 
ring clumps, with an average observed Tb of ^5 K and 
LVG predicted T b of ~20 K, then f b is -0.25 and 9 S 
~44 pc. The size of the broad line ring clumps is a fac- 
tor of 2 higher than the narrow line ring clumps. Given 
the estimated low volume number density of the broad 
line clumps in Sect. l3.4~Tl the inconsistency between high 
Sh 2 and low number density mentioned is possibly due 
to the large line-of-sight path, or larger scale height of 
the broad line clumps. If the assumption of a spherical 
shape of the clumps hold, the larger size estimated above 
seems to be consistent with this scenario. However, con- 
sidering the higher opacity of the broad line clumps as 
suggested by the high £h 2 ; h is likely that we are trac- 
ing the diffuse and cold gas at the surface of GMA rather 
than the dense and warm gas. On the other hand, the 
narrow line clumps show the opposite trends to trace the 
relatively warm and dense gas. 

4.1.3. Azimuthal variation of the line widths and Eh 2 

In Figure [51 we found the observed (and intrinsic) line 
widths, and £h 2 , show variations along the azimuthal 
direction in the ring. There are local maxima of line 
widths and Sh 2 at the position of the twin-peak at the 
orbit crowding regions. Since the brightness tempera- 
ture does not vary dramatically between the emission 
peaks along the ring, the deduced peaks in H 2 column 
density can be directly attributed to the increased line 
widths. Beam smearing is possibly important as mul- 
tiple streams may converge within a synthesized beam. 



However, as shown in Fig lllb . enhanced line widths oc- 
cur even over extended regions. The intrinsic velocity 
dispersion of both the narrow and the broad line clumps 
cannot be thermal, as the implied temperature would 
be 2000-10000 K. In the Galactic GMCs, non-thermal 
line broadening has been attributed to turbulence from 
unknown mechanisms. In the case of NGC 1097, turbu- 
lence could have been generated by shocks of the orbit 
crowding region. Shocked gas as indicated by H 2 S(l-0) 
emissi on has been reported to ward the twin-peak of NGC 
1097 (Kotilain en et al.l 120001) . The narrow line clumps 
further along the ring may be due to subsequent dissi- 
pation of the kinetic energy from the shock wave. This 
kind of large velocity widths up to ~100 km s _1 were 
also observed in NGC 6946. With resolutions down to 
CMC-scale of 10 pc, multiple compone nts in the spectra 
were s een at the nuclear twin-peak by ISchinnerer et all 
(|2007[) . However, the mechanism to cause the broadened 
line widths is not clear at this scale. 

As noted earlier, the "intrinsic" line widths we derived 
contain not only the random dispersion but also the ve- 
locity gradients due to non-circular motions generated 
by shock fronts. However, our observations do not have 
sufficient angular resolution to resolve the locations and 
magnitudes of the discre te velocity jumps to be expected 
across the shock front ([Draine fc McKed[l993l) . A de- 
tailed hydrodynamical model is also needed to quantita- 
tively predict the magnitude of the non-circular velocity 
gradient in the ring. 

4.2. Gravitational stabilities of the GMAs in the 
starburst ring 

How do the GMAs form in the starburst ring? In the 
ring of NGC 1097, we c onsider gravit ational collapse due 
to Toomre instability (|Toomrel |i964'). We estimate the 
Toomre Q parameter (S C rit/^gas) to see if the rotation 
of the ring is able to stabilize against the fragmentation 
into clumps. Here the Toomre critical density can be 
expressed as: 



a 



V 

k = 1.414 | — 
R 



3.36G/ ' 



RdV\ 
1+ VdR) 



0.5 



(11) 



(12) 



where the constant a is unity, G is the gravitational con- 
stant, k is the epicycle frequency, V is rotational veloc- 
ity, and R is the radius of the ring. If Eh 2 exceeds S C rit, 
then the gas will be gravitationally unstable and collapse. 
We approximate the velocity gradient dV/dR as close to 
zero in Figure Q21 because of the flat rotation curve of 
the galaxy at the position of the starburst ring. There- 
fore K ~ 1.414V7i?. V is ~ 338 km s" 1 in the galactic 
plane assuming an inclination of ~42°. We find the ratio 
of S cr it/SH 2 is less than unity in the ring (~0.6). Since 
S gas consists of H2, HI, and metals, this ratio of 0.6 is 
an upper limit. However, the HI is often absent in cen- 
ters of galaxie s, and NGC 1097 also shows a HI hole in 
the central 1' (jHigdon fc Wa llin 2003). Hence the HI gas 
does not have an important contribution in the nuclear 
region. A ratio less than unity suggests that the ring is 
unstable and will fragment into clumps. 
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As for whether the GMA itself is gravitational bound 
or not, in Sect. 13 .2.21 we found that M v i r /Af gas is around 
unity in the narrow line clumps and larger than unity 
by a factor of more than 2 in the broad line clumps. 
This seems to suggest that the broad line clumps are 
not virialized, probably because of the larger turbulence. 
However, there are several factors to reduce the ratio. 
First, since we do not subtract the non-circular mo- 
tion in the broad line clumps, the observed M vir /M gas 
should be an upper limit. We can estimate how large 
the non-circular motion is if we assume the M V i r /A/ gas 
of the broad line clumps is unity. The magnitude of the 
non-circular motion is from 50 to 100 k m s" 1 for the 
broad l ine clumps under this assumption. lAthanassoulal 
(|1992al) showed that there is a correlation between the 
bar axial ratio and the velocity gradient (jump) across 
the s hock wave. NGC 1097 has a b ar axial ratio of 
~2.6 (Mencnd ez-Delmestre et all 120071 ). which indicates 
the maximum velocit y jump is ~ 70 km s" 1 (Figure 12; 
lAthanassoulal Il992bl ) . However, this number is sup- 
posed to be an upper limit since it was measured at the 
strongest strength of the shock, where the strength of 
the shock is a function of position relative to the nu- 
cleus. The shock strength seems to be weaker at the 
intersection of the circumnuclear ring than for the outer 
straight dust lanes, as suggested in the model. Hence 
we expect the velocity gradient caused by shock front is 
smaller than ~ 70 km s _1 , based on this correlation. In 
this case, the broad line clumps still have 20% larger ve- 
locity dispersion than narrow line clumps. Second, the 
size of the clumps will also be a possible factor to reduce 
the M v ; r /M gas to unity. In Sect. 14. Ol we estimate that 
the filling factor of the broad line clumps is —0.25, and 
therefore the intrinsic radius will be smaller by a factor 
of 0.25 1 / 2 -0.5. Third, in Sect. ET2~2l we point out the 
molecular gas mass derived from 12 CO (J = 2-1) might 
be underestimated at least by a factor of ~2. These fac- 
tors also can lower the M v i r / M gas to roughly unity in 
most of the broad line clumps and hence the GMA could 
also be gravitationally bound in the broad line clumps. 

4.3. Star formation in the ring 
4.3.1. Extinction 

The distribution of the massive star clusters is uni- 
form in the ring instead of highly clustering in certain 
clumps in Figure [TJ The star clusters do not coincide 
with most of the CO peaks. This could be due to several 
reasons such as extinction or the physical nature of the 
star clusters. Of course, since we are comparing the scale 
of star clusters (2 pc) with that of GMAs (100 pc), we 
are not able to conclude the physical correlation by their 
spatial distribution. The Paa is more transparent than 
the commonly used Ha through dust extinction. How- 
ever, the foreground extinction (E(B-V) — 1.3) for the 
Paa clusters corre sponds to a hydrogen c olumn density 
of 6.4x 10 cm" " 2 (jDiplas fc Savagdfl99l averaged over 
one synthesized beam. This is much less than the average 
H2 column density for the molecular clumps in the ring, 
which is — 8.7xl0 22 cm" 2 . Therefore, this suggests that 
the detected Paa clusters might be located on the sur- 
face of the clouds instead of being embedded inside the 
clumps, or else located away from the clouds. However, 
it can not be ruled out that there are deeply embedded 



stellar clusters in the CO peaks in this scenario. 

There is also a deficiency of star clusters in the broad 
line clumps associated with the dust lanes. This again 
could be due to extinction although we found it is similar 
in the star forming ring and the dust lanes. However, it 
is interesting that there is a spatial offset between the 
Spitzer 24 /im peaks and the CO peaks, and the FIR 
emission in the dust lane is intrinsical ly faint based on 
Hersc hel PACS 70 and 100 fim maps (jSandstrom et al.1 
2010). The long wavelength IR results are less affected 
by extinction, and suggest a lack of newly formed star 
clusters in the dust lanes. However, a higher resolution 
for the FIR observation is needed to confirm the star 
formation activities in the broad line clumps. 

4.3.2. Molecular gas and star formation 

In NGC 1097, the mechanisms of the intensive star for- 
mation in the ring are still uncertain. It could be induced 
by the grayitatio nal collapse in the ring stochastically 
(Elmegreen 1994). The other possible mechanism is that 
the stars form in the downstream of the dust lane at the 
conjunction of the ring, an d the star clusters continue to 
orbit along the ring (e.g.. iBoker et all 120 08). The ma- 
jor difference is that the latter scenario will have an age 
gradient for the star clusters along the ring while it is 
randomized in the previous case. Several papers have 
discussed these mechanisms in the galaxies that have 
star forming rings and there is no clear answer so far 
(e g iMazzuca et al.ll2008HBoker et al.ll2008t iButa et al.1 
2000). ISandstrom et al.l 1)20101 ) tested the above pictures 
by examining if there is an azimuthal gradient of dust 
temperature in the ring of NGC 1097, assuming that the 
younger population of massive star clusters will heat the 
dust to higher temperatures than older clusters. There 
seems to be no gradient, though it is difficult to conclude 
since a few rounds of galactic rotation might smooth out 
the age gradient. We do not aim to solve the above ques- 
tion in this paper since the most direct way is to measure 
the age of the star clusters, and this needs detailed mod- 
eling. It is interesting to compare the properties of the 
molecular clumps with star clusters since the molecular 
clouds are the parent site of star formation. The one rel- 
evant result here is that the line widths are narrower fur- 
ther away from where the molecular arms join the ring. 
This could be related to the dissipation of turbulence 
which may allow cloud collapse to proceed. 

In Figure[9jb), we found that Ssfr, compared with ve- 
locity dispersion (Figure [8]), has no significant azimuthal 
correlation. Furthermore, R32 shows a similar trend to 
Ssfr as a function of azimuthal direction (Figure [S^d)), 
where R32 and Ssfr has a correlation in Figure [9^c). 
This suggests that Ssfr an( i are physically related, 
but might not be associated with the large scale dynamics 
in this galaxy. Although Ssfr seems to be suppressed in 
the broad line clumps. Nevertheless, in Figure [§Ib)&(d) 
we consider that the standard deviations of the measured 
Ssfr do not significantly diverse among global varia- 
tions. The standard deviations are 0.7 M Q yr" 1 kpc" 2 
and 0.4 M yr" 1 kpc" 2 of the northeast and southwest 
of the ring, respectively. The mean values of the Ssfr 
are 3.3 M© yr _1 kpc -2 and 1.4 M Q yr" 1 kpc" 2 of the 
northeast and southwest of the ring, respectively. With 
the limited points, our results suggest that the star for- 
mation activities are randomly generated on the local 
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scale instead of a systematic distribution. 

In Figure [13l we show how the R 32 ratio varies for 
different densities and kinetic temperatures. It shows 
that when R32 varies from 0.3 to 0.9, the required number 
density of molecular gas changes from 10 2 cm -3 to 5x 10 3 
cm~ 3 . The R32 seems to be dependent on the density 
more than the temperature when it is below unity. Hence 
the variation of R32 indicates different density among 
clumps. It is interesting to note that there is a correlation 
between R 32 and SFR, and some clumps (Nl, N2, N7) 
which have higher SF R values are spatia lly close to the 
HCN(J = 1-0) peaks (|Kohno et al.ll200l . Higher values 
of R32 will select denser gas associated with higher SFR, 
as has been shown in the large scale observations of HCN 
and FIR correlation (iGao. Solomon, fc Philip) I2004D . In 
the smaller GMC-scale. lLadal (|1992| ) also showed that 
the efficiency of star formation is higher in the dense 
core rather than in the diffuse gas. The ratio map of R32 
can be useful to determine the location of star formation. 
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5. SUMMARY 

1. We show the multi-J CO line maps of NGC 1097 to- 
ward the 1 kpc circumnuclear region. The molecu- 
lar ring is resolved into individual GMAs in the star 
forming ring, the dust lanes, and at the twin-peak 
structures. For the first time the molecular concen- 
tration at the twin-peak is resolved in to two pop- 
ulations of GMAs in terms of velocity dispersion 
and physical conditions. The clumps in the star- 
burst ring have narrower velocity dispersion while 
the line widths are broader in the dust lanes, and 
for some clumps located in the twin-peak. The 
physical and kinematic properties are different for 
these clumps. The narrow line clumps have higher 
temperature (>250 K) and density ((4.5±3.5) x 10 3 
cm -3 ) in contrast to the broad line clumps (T = 
45±15 K; N H , 2 = (8.5 ± 1.5) x 10 2 cm" 3 based on 
the LVG analysis. 



2. The Toomre-Q factor is smaller than unity in the 
molecular ring suggesting that the GMAs could 
form via gravitational instability in the ring, where 
the Sh 2 of the clumps is large enough to over- 
come the critical density. The narrow line clumps 
are gravitationally bound as shown by the values 
of Mvir/Mgas which are nearly unity. Although 
Mvir/Mgas is larger than unity in the broad line 
clumps, by accounting for non-circular motions, 
smaller intrinsic source sizes, and the underestima- 
tion of molecular gas mass, we can lower M v j r / M gas 
to unity. Therefore both systems are likely to be 
gravitationally bound. 



3. The SFR is correlated to R32, suggesting that the 
star formation activities and the physical condi- 
tions of the molecular gas are associated with each 
other. In contrast to the velocity dispersion and 
the Sh 2 , the SFR and R32 are not correlated with 
the large scale dynamics. This suggests that the 
visible star formation activities remain a localized 
phenomenon. The SFR is lower in the broad line 
clumps than in the narrow line clumps, which may 
be intrinsicly suppressed in the dust lanes. 



CO(J = 2-1) Imaging of Seyfert 1 Galaxy NGC 1097 



Table 1 

SMA observation parameters 





230 GHz 


230 GHz 


230 GHz 


345 GHz 


Parameters 


Compact-N 


Extended 


Very extended 


Compact 


Date 


2004-07-23, 2004-10-01 


2005-09-25 


2005-11-07 


2006-09-05 


Phase center (J2000.0): 










R.A. 




<*2000 — 






Dccl. 




^2000 — 


-30°16'28 "897 




Primary beams 




52" 




36" 


No. of antennas 


8, 8 


6 


7 


7 


Project baseline range (kA) 


5 - 74, 10 - 84 


23 - 121 


12 - 390 


7.2 - 80 


Bandwidth (GHz) 






1.989 




Spectral resolution (MHz) 


0.8125, 3.25 


0.8125 


0.8125 


0.8125 


Central frequency, LSB/USB (GHz) 




219/228 




334/344 


T225 b 


0.15, 0.3 


0.06 


0.1 


0.06 


T sys ,DSB (K) 


200, 350 


110 


180 


300 


Bandpass calibrators 


Uranus, J0423 - 013 


3C454.3 


3C454.3, 3C111 


Uranus, Neptune 


Absolute flux calibrators 3 


Uranus (36.8, 34.2), 


Uranus (37.2), 


Uranus (34.9), 


Neptune (21.1) 




J0423 - 013 (2.8, 2.5) 


3C454.3 (21.3) 


3C454.3 (21.3) 


Gain calibrators 


J0132 - 169 


J0132 - 169 


J0132 - 169 


J0132 - 169, J0423 - 013 



a The numbers in the parenthesis arc the absolute flux in Jy. 
b T225 is the optical depth measured in 225 GHz. 
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Table 2 

Physical parameters of the peaks of the molecular clouds 



ID 


SR. A. 

I'll 

V ) 


<5Dccl. 
(2) 
V ) 


(3) 

(Jy bccim km s ) 


SV u 

(4) 

1 

(km s ) 


SV- t 

u * mt 

(5) 

1 

(km s ) 


iv H 2 
(6) 


(7) 

l lvi O P c ; 


Ma 
(8) 

(iv ivi0; 


± b 


Nl 


9.0 


1.6 


17.2 


54±3 


42±4 


7.9 


1270 


10.7 


4.1 


N2 


8.6 


4.8 


27.9 


67±3 


61±4 


12.8 


2070 


17.3 


5.5 


N3 


3.6 


8.4 


38.3 


69±1 


64±1 


17.6 


2840 


23.8 


7.5 


N4 


0.6 


8.2 


18.5 


62±4 


56±4 


8.5 


1370 


11.5 


5.8 


N5 


-4.8 


7.8 


21.3 


38±1 


31±1 


9.8 


1580 


13.2 


6.1 


N6 


-6.4 


5.6 


19.6 


41±2 


37±2 


9.0 


1460 


12.2 


4.1 


N7 


-8.6 


-2.6 


40.3 


61±2 


52±3 


18.5 


2990 


25.0 


8.0 


N8 


-5.4 


-8.4 


28.4 


57±2 


51±2 


13.1 


2110 


17.6 


5.8 


N9 


-2.8 


-9.2 


24.6 


62±3 


57±3 


11.3 


1820 


15.2 


4.1 


N10 


4.8 


-7.8 


19.7 


43±2 


37±2 


9.1 


1460 


12.2 


6.1 


Nil 


7.8 


-6.6 


27.7 


52±4 


49±4 


12.7 


2050 


17.2 


5.8 


Bl 


10.2 


3.6 


27.8 


113±10 


109±11 


12.7 


2060 


17.2 


2.4 


B2 


-7.4 


-6.8 


34.8 


99±5 


96±5 


16.0 


2580 


21.6 


4.6 


B3 


-5.4 


-6.6 


31.3 


94±8 


90±8 


14.4 


2320 


19.4 


4.1 


Dl 


16.8 


-1.4 


18.9 


100±10 


97±10 


8.7 


1400 


11.7 


1.8 


D2 


12.8 


1.6 


21.1 


86±17 


82±18 


9.7 


1560 


13.1 


2.0 


D3 


-17.8 


2.6 


16.2 


81±7 


78 ± 8 


7.5 


1200 


10.1 


1.8 


D4 


-13.0 


-1.0 


24.4 


84±5 


80±6 


11.2 


1810 


15.1 


3.5 


D5 


-10.0 


-4.0 


38.6 


118±4 


115±4 


17.7 


2860 


24.0 


4.6 


Nu 








50.4 


52±5 




16.1 


3740 


23.2 




Nu 










57±5 












Nu 










186±20 













Note. — We define the peaks based on their location and their velocity dispersions. The clumps in the dust lanes (molecular 
spiral arms) arc named Dl, D5. The clumps in the ring are further designated by their velocity dispersion being broader or 
narrower than 30 km s — 1 , and named respectively as Bl, B3, and Nl, Nil. Nu is the ID of the nucleus. (1) R.A. offsets 
from the phase center. (2) Dec. offsets from the phase center. (3) Integrated CO(J — 2 — 1) intensity. The uncertainty is 2.2 
Jy beam -1 km s — 1 . (4) Fitted FWHM for the observed line width. The nucleus has multiplc-gaussians profile, and we list the 
fitted line widths with 3 gaussians. Note that the IqO: ^h 2 j an d Sh 2 of the nucleus is the sum value of the three components. 
(5) FWHM line width of intrinsic velocity dispersion. (6) H2 column density with the uncertainty of l.lxlO 22 cm -2 . (7) Eh 2 
with the uncertainty of 170 Mq pc -2 . (8) Mass of molecular H2 within the synthesized beam (l'.'5 X l^O). The uncertainty is 
1.4 Xl0 6 Mq. (9) Peak brightness temperature. 



CO(J = 2-1) Imaging of Seyfcrt 1 Galaxy NGC 1097 



Table 3 

Physical parameters of the GMAs 



ID 


Diameter 


M H2 
(2) 


-Mgas 


M ■ 


M vir /M gas 












(5) 




( ) 


(10 6 M ) 










Nl 


2.9 


32.4 ± 


5.3 


44.1 ± 7.2 


44.9 ± 10.2 


1.0 


± 0.3 


N2 


3.7 


71.2 ± 


8.7 


96.9 ± 11.8 


118.0 ± 16.3 


1.2 


± 0.2 


N3 


3.5 


73.6 ± 


7.9 


100.1 ± 10.8 


125.4 ± 5.9 


1.3 


± 0.1 


N4 


2.2 


17.7 ± 


3.0 


24.0 ±4.1 


58.6 ± 9.6 


2.4 


± 0.6 


N5 


2.9 


37.5 ± 


5.5 


51.0 ± 7.4 


25.1 ± 2.7 


0.5 


± 0.1 


N6 


2.9 


33.8 ± 


5.6 


46.0 ± 7.6 


34.6 ± 5.2 


0.8 


± 0.2 


N7 


3.4 


79.1 ± 


7.6 


107.5 ± 10.3 


79.3 ± 9.9 


0.7 


± 0.1 


N8 


2.5 


51.5 ± 


4.2 


70.0 ± 5.7 


57.5 ± 5.6 


0.8 


± 0.1 


N9 


3.7 


94.5 ± 


8.8 


128.5 ± 12.0 


104.2 ± 14.1 


0.8 


± 0.1 


N10 


2.8 


32.7 ± 


5.2 


44.5 ± 7.1 


34.3 ± 4.1 


0.8 


± 0.2 


Nil 


4.1 


133.1 ± 


11.1 


181.0 ± 15.0 


85.8 ± 15.1 


0.5 


± 0.1 


Bl 


3.1 


48.6 ± 


6.3 


66.1 ± 8.6 


316.6 ± 71.0 


4.8 


± 1.2 


B2 


3.3 


80.0 ± 


7.2 


108.8 ± 9.9 


267.4 ± 29.8 


2.5 


± 0.4 


B3 


3.0 


77.0 ± 


5.7 


104.7 ± 7.7 


208.8 ± 43.4 


2.0 


± 0.4 


Dl 


3.3 


56.5 ± 


6.9 


76.9 ± 9.4 


269.1 ± 62.8 


3.5 


± 0.9 


D2 


3.5 


58.6 ± 


8.2 


79.7 ± 11.1 


207.7 ± 102.0 


2.6 


± 1.3 


D3 


2.9 


30.1 ± 


5.5 


41.0 ± 7.4 


152.3 ± 34.2 


3.7 


± 1.1 


D4 


5.4 


117.7 ± 


18.9 


160.0 ± 25.7 


303.0 ± 47.2 


1.9 


± 0.4 


D5 


4.1 


112.6 ± 


11.0 


153.1 ± 15.0 


474.8 ± 36.8 


3.1 


± 0.4 



Note. — (1) Diameter of the clumps. (2) H2 mass integrated over the diameter 
of the clumps. (3) Gas mass integrated over the diameter of the clumps. The M gas 
is the Mh 2 corrected with the He fraction of 1.36 in Sect.fefscct-mass. (4) Virial 
mass of the clumps. (5) Ratio of the virial mass to the M gas . We have not corrected 
for the beam convolution effect as the derived diameters arc sufficiently large. We 
think the assumption of Gaussian shapes may be the greater source of error. 
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Table 4 

Star formation properties 



ID 




B(B-V) 


A 


•Paa 




A v 


SSPR, 




(!) . 


(2) 




(3) 




(4) 


(5) 






(10 14 erg s 1 cm 2 ) 


(mag) 


(mag) 


(mag) 


/TV J I 

(Mq yr 


kpc J 


Nl 


6.03 ± 0.55 


1.41 ± 0.01 


0.64 


± 0.30 


4.36 


± 0.04 


2.74 ± 


0.25 


N2 


7.01 ± 0.39 


1.06 ± 0.01 


0.48 


± 0.23 


3.30 


± 0.04 


3.18 ± 


0.18 


N3 


4.00 ± 0.22 


1.02 ± 0.02 


0.46 


± 0.22 


3.16 


± 0.06 


1.81 ± 


0.10 


N4 


8.95 ± 0.62 


1.21 ± 0.01 


0.55 


± 0.26 


3.75 


± 0.03 


4.06 ± 


0.28 


N5 


7.52 ± 0.38 


1.00 ± 0.01 


0.46 


± 0.22 


3.11 


± 0.03 


3.41 ± 


0.17 


N6 


7.94 ± 0.37 


0.97 ± 0.01 


0.44 


± 0.21 


3.00 


± 0.03 


3.60 ± 


0.17 


N7 


8.97 ± 0.75 


1.36 ± 0.01 


0.62 


± 0.29 


4.20 


± 0.03 


4.07 ± 


0.34 


N8 


4.21 ± 0.41 


1.45 ± 0.02 


0.66 


± 0.31 


4.50 


± 0.06 


1.91 ± 


0.19 


N9 


3.72 ± 0.35 


1.41 ± 0.02 


0.64 


± 0.30 


4.36 


± 0.07 


1.69 ± 


0.16 


N10 


2.90 ± 0.30 


1.47 ± 0.03 


0.67 


± 0.32 


4.56 


± 0.09 


1.32 ± 


0.14 


Nil 


1.92 ± 0.22 


1.49 ± 0.05 


0.68 


± 0.32 


4.62 


± 0.14 


0.87 ± 


0.10 


Bl 


2.55 ± 0.19 


1.17 ± 0.03 


0.53 


± 0.25 


3.62 


± 0.10 


1.16 ± 


0.09 


B2 


2.60 ± 0.24 


1.35 ± 0.03 


0.61 


± 0.29 


4.18 


± 0.10 


1.18 ± 


0.11 


B3 


1.29 ± 0.17 


1.44 ± 0.07 


0.65 


± 0.31 


4.45 


± 0.21 


0.58 ± 


0.08 


Dl 


<0.12 
















D2 


<0.12 
















D3 


<0.12 
















D4 


<0.12 
















D5 


1.56 ± 0.16 


1.31 ± 0.05 


0.60 


± 0.28 


4.07 


± 0.17 


0.71 ± 


0.07 



Note. — (1) Pact flux corrected by the extinction measured in the CO clumps. The upper 
limit of the clumps (Dl, .... D4) is 1.21 xl0~ 15 erg s -1 cm -2 . (2) Color excess. (3) Extinction at 
wavelength of Pact in the unit of magnitude. (4) Extinction at V-band in the unit of magnitude. (5) 
Surface density of SFR. The upper limit is 0.05 Mq yr -1 kpc -2 . 



CO(J = 2-1) Imaging of Seyfert 1 Galaxy NGC 1097 



Table 5 

Intensity ratios of the CO clumps 



ID 


R32 


Rl3 




(l) 


(2) 


Nl 


1.18 ± 0.07 (0.78 ± 0.18) 


10.52 ± 7.05 


N2 


0.89 ± 0.03 (1.02 ± 0.14) 


9.17 ± 2.92 


N3 


1.22 ± 0.04 


8.12 ± 2.35 


N4 


1.41 ± 0.09 


- 


N5 


1.68 ± 0.09 


12.44 ± 9.74 


N6 


1.19 ± 0.07 


12.47 ± 9.50 


N7 


0.92 ± 0.02 (0.99 ± 0.06) 


8.52 ± 1.71 


N8 


0.62 ± 0.02 (0.50 ± 0.06) 


10.24 ± 3.17 


N9 


0.66 ± 0.03 


10.88 ± 5.14 


N10 


0.85 ± 0.05 


7.18 ± 2.73 


Nil 


0.78 ± 0.03 


13.04 ± 6.45 


Bl 


0.79 ± 0.03 




B2 


0.71 ± 0.02 


6.51 ± 1.10 


B3 


0.72 ± 0.02 


11.49 ± 3.74 


Dl 


0.76 ± 0.08 




D2 


0.86 ± 0.05 




D3 


1.19 ± 0.12 




D4 


0.72 ± 0.03 




D5 


0.72 ± 0.02 (0.39 ± 0.03) 


10.65 ± 2.59 (8.38 ± 3.41) 


Nu 


0.93 ± 0.02 


23.37 ± 9.79 



Note. — (1) 12 CO(J = 3-2)/(J = 2-1) intensity ratios derived 
from CO brightness temperature. Numbers within parenthesis arc 
beam smearing corrected ratios. (2) 12 CO(J = 2-l)/ 13 CO(J = 2- 
1) intensity ratios derived from CO brightness temperature. Num- 
bers within parenthesis are beam smearing corrected ratios. All the 
quantities are measured with a beam size of 3''25x2''55. 



Table 6 

Dynamical parameters fitted by GAL 



Parameters Value 

R.A. 02 h 46 m 18?95 

Decl. -30°16'29^'13 

Position angle 133?0±0?1 

Inclination 41?7±0?6 

Systemetic velocity (km s _1 ; V sys ) 1249.0±0.5 

Vmax (km s- 1 ) 387.6±4.3 

Kmax 9'.'5±0'.'1 

Note. — The 6-cm peak of the nucleus is R.A. — 
02 h 46 m 18?96, Decl. = -30° 16'28"897. 
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Figure 1. Top left image is the 12 CO(J = 2-1) integrated intensity map. The contour levels are 2, 3, 5, 20, 25, and 30cr (1<7 = 2.3 
Jy km s" 1 beam- 1 ). The synthesized beam is l'.'5xl'.'0 (PA = 8?1). The emission at (10", -15") are sidelobes. Top right image is the 
12 CO(J = 3-2) integrated intensity map. The contour levels are 5, 7, 9, 10, 15,..., 60, 80, and IOOct (ler = 3.0 Jy km s" 1 beam" 1 ). The 
synthesized beam is 3"5x2 / /l (PA = -4? 4). Bottom left image is the 13 CO(J = 2-1) integrated intensity map. The contour levels are 2, 3, 
and 4cr (la = 2.1 Jy km s - 1 beam" 1 ). The synthesized beam is l'.'8xl'.'4 (PA = 19°). Bottom right image is the 12 CO(J = 2-1) peak 
brightness temperature map. The contours are 2, 3, 4, 5, 6, and 7 K. All of the maps are overlaid with the positions of Pact star clusters 
(squares), 6 cm radio continuum sources (crosses), and V-band (< 13 mag) star clusters (circles). The central cross in each map is the 
position of the 6-cm nucleus (Hummel et al. 1987), which is assumed to be the active nucleus. The beam size is shown in the lower right 
corner of each map. 
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Figure 2. The 12 CO(J = 2-1) channel maps are overlaid on the archival HST I-band (F814W) image with the corrected astrometry. 
The contour levels are —2, 2, 4, 8, 16, and 32cr, where la = 20 mjy beam -1 (306 mK) in 40 km s _1 resolution. The velocity (km s — 1 ) with 
respect to the systemic velocity of 1254 km s _1 (Kohno ct al. 2003) is labeled in the top left corner of each map. The beam size (l'.'5xl'.'0, 
PA = 8?1) is shown in the lower right corner of each map with solid ellipse. The dirty beam is at the bottom right panel with a contours 
level of -100, -50, -10, -5, 5, 10, 50, 100% of the peak. 
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Figure 3. The 12 CO(J = 3-2) channel maps are overlaid on the archival HST I-band (F814W) image with the corrected astrometry. 
The contour levels are -3, 3, 5, 10, 20, 40, 60, and 80cr, where la = 17.5 mjy beam - (23 mK) in 40 km s -1 resolution. The velocity (km 
s" 1 ) with respect to the systemic velocity of 1254 km s -1 is labeled in the top left corner of each map. The beam size (3''5x2'.'l, PA = 
—4? 4) is shown in the lower right corner with solid ellipse. The dirty beam is at the bottom right panel with a contours level of —100, —50, 
-10, -5, 5, 10, 50, 100% of the peak. 
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Figure 4. Top image is the 12 CO(J = 2-1) integrated map (contours) overlaid on the archival HST I-band (Filter F814 W) image (grey 
scale). Astrometry of the HST image was corrected using background stars with known positions. The contour levels for the 12 CO(J = 
2—1) are 2, 3, 5, 20, 25, and 30 a (1 a = 2.3 Jy km s — 1 beam -1 ). The IDs for the individual peaks of clumps are marked. The CO 
synthesized beam (l'.'5x l'.'O) is shown in the lower right corner. Bottom image is the HST NICMOS Pact line image (color) overlaid on 
the 12 CO(J = 2-1) contour. The contour levels are the same as in upper image. 
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Figure 5. The high resolution (l'.'5x l'.'O) 12 CO(J = 2—1) spectra of individual clumps measured at the 12 CO(J = 2-1) peak position 
within one beam. The velocity is relative to the systemic velocity of 1254 km s — 1 and spectral resolution is 10 km s — 1 . The IDs of the 
clumps are labeled in each panel. 
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Figure 6. Left: We show the 12 CO(J = 2-1) spectra, where the solid line is the JCMT data (Petitpas et al. 2003) and the dotted line 
is our SMA data. The intensity scale is the main beam temperature at 21" resolution. We only took 6 chunks to make the SMA maps, 
so the velocity range is smaller than that of JCMT. Right: The 13 CO(J = 2-1) spectra, where the solid line is the JCMT data, and the 
dotted line is our SMA data. The beam size of the two data are matched to 21". 
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Figure 7. (a) The number histogram of the total H2 mass of the clumps in the ring. Horizontal and vertical axis are the gas mass and 
number, respectively. The gas mass is in units of 10 6 Mq. The negative horizontal axis is to show the plot clearly, (b) The correlation 
between total H2 mass and FWHM intrinsic line width of the narrow line ring clumps (circles) , broad line ring clumps (squares) , and dust 
lane clumps (triangles). The H2 mass is in units of 10 6 Mq. The uncertainties of ±lcr are overlaid on the symbols with vertical/horizontal 
bars. 
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Figure 8. (a) Sh 2 measured at the position of the intensity peak in units of Mq pc -2 (see Table [2} as a function of azimuthal angle. 
East direction corresponds to 0°, and increase is in clockwise direction. The dashed lines mark position angle from 0° to 45° and from 
180° to 225°, which roughly correspond to the position of the orbit crowding regions. The meaning of the symbols are the same as in 
Figure [7p. (b) The FWHM intrinsic line width of clumps as a function of azimuthal angle. The dashed lines are the same as in (a). The 
uncertainties of ±1<t are overlaid on the symbols with vertical/horizontal bars. 
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Figure 9. (a) Surface SFR density (Mq yr — 1 kpc -2 ) is shown as a function of Ejj 2 (Mq pc — 2 ). The symbols are the same as in Figurel8l 
The meaning of the symbols are the same as in Figure [7p. (b) Surface SFR density is shown as a function of azimuthal angle. The dashed 
lines mark position angle from 0° to 45°and from 180° to 225°. (c) Surface SFR density (Mq yr _1 kpc -2 ) is shown as a correlation of 
R32- (d) R32 is shown as a function of azimuthal angle. The uncertainties of ±ltr are overlaid on the symbols with vertical/horizontal bars. 
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Figure 10. The correlation of Sjj 2 and surface SFR density of individual clumps of NGC 1097 are overlaid on the data used in Kennicutt 
(1998). Their data for normal galaxies are represented as plus signs, and infrared-selected circumnuclear starburst galaxies as asterisks, 
while their NGC 1097 data point is marked as triangle. Our spatially resolved clumps of the circumnuclear starburst ring of NGC 1097 are 
marked as crosses, and the average value of the clumps is marked as a square. 
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Figure 11. (a) The intensity weighted mean velocity map (MOM1) of 12 CO(J = 2-1) line with respect to the systematic velocity (1254 km 
s _1 ), solid and dashed lines represent the redshifted and blueshifted velocity respectively. The first negative contour (close to the central 
cross) is km s — 1 , and the contour spacing is in 25 km s — 1 resolution, (b) The intensity weighted velocity dispersion map (MOM2). The 
contour interval is 10 km s _1 , note the values are not FWHM line width but the square root of the dispersion relative to the mean velocity. 
Therefore the number is lower than the FWHM we derived by line fitting. 
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(a) projected view (b) intrinsic view 




Figure 14. The cartoon sketch of the gas morphology in the circumnuclcar region of NGC 1097. The red circle represents the starburst 
ring, where the narrow line clumps are located. The black curves are the dust lane associated with shock wave, where the broad line 
clumps are located. The blue line is the major axis of the large scale stellar bar. We show: (a) the projected view of the morphology of 
the starburst ring/dust lane associated with our observation. It shows a nearly circular starburst ring, (b) the intrinsic shape of the ring, 
which is expected to be an ellipse. 



